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Kinematic distributions in the decays of the newly discovered resonance to four leptons can provide 
a direct measurement of the tensor structure of the particle's couplings to gauge bosons. Even if the 
particle is shown to be a parity even scalar, measuring this tensor structure is a necessary step in 
determining if this particle is responsible for giving mass to the Z. We consider a Standard Model like 
coupling as well as coupling via a dimension five operator to either ZZ or Zj. We show that using 
full kinematic information from each event allows discrimination between renormalizable and higher 
dimensional coupling to ZZ at the 95% confidence level with 0(50) signal events, and coupling to 
Z'y can be distinguished with as few as 20 signal events. This shows that these measurements can 
be useful even with this year's LHC data. 



I. INTRODUCTION 

The discovery of a new particle by the CMS [U 

and 

ATLAS experiments is a great triumph for particle 
physics, but it also leads to a host of new questions about 
the nature of the new state. The question of whether this 
particle gives mass the the W and Z is of paramount 
importance. Since the discovery, the theoretical commu- 
nity has begun to weigh in on this question with many 
different analyses [3-26]. These analyses show that the 
new state is broadly consistent with the Standard Model 
Higgs, but all of these analyses use only the informa- 
tion obtained from cross sections and branching ratios. 
While this form of analysis is very useful, it is ultimately 
a model dependent test of the properties of the new state. 

Here we attempt to delve into the signal events them- 
selves and see what information can be learned about the 
new particle from the kinematic distributions of the final 
state particles. Specifically, the channel where the new 
particle decays to four leptons via intermediate gauge 
bosons [27L 12a ] contains a tremendous amount of infor- 
mation about the new resonance. It has been shown that 
this channel can be used to distinguish the parity and 
spin (29l - [37j of a new resonance, with [35| and [36| do- 
ing an exhaustive comparison of many different spin and 
parity possibilities. It has also been shown that kine- 
matic methods can be used to distinguish signal from 
background in this channel (38|. 

Here we take the hypothesis that the new particle is 
a parity even scalar and try to see if this channel can 
be used to directly measure the tensor structure of the 
coupling of this particle to the four lepton final state. If 
we denote the new scalar by <f>, it can have the following 
couplings to ZZ 



a h m%> 



<Z„Z» 



>Z^Z 



(1) 



Here 



where Z M is the Z field while Z^ — d^Z^ 
v = 246 GeV is the Standard Model (SM) Higgs vev 
which is chosen to normalize the operators, and the ... 



is for operators of dimension higher than five. If 4> is 
the Standard Model Higgs, then ah — i, and the other 
coupling is loop induced and small. 

As we are trying to determine whether this new parti- 
cle is the SM Higgs, we must consider other possibilities. 
If <j) does not give mass to the Z, then its linear coupling 
to gauge bosons can proceed via the field strength tensor, 
Z^ v as in the operator a s in Eq. ([1} . There are many such 
models in the literature, see for example j3f| [3^, H3] and 
references therein. The a s operator is generically loop 
induced and its coefficient is model dependent. We see 
this sort of operator even in the Standard Model Higgs' 
coupling to 77 and Zj: 



(2) 



where we continue to use <j> to denote our scalar, and F^ 
is the field strength tensor for the photon. In the SM, 
a 7 and az-y are induced by loops with top and W giving 
the largest contributions. If cf> is not the Higgs, then a 
plausible alternative is that it decays to four leptons via 
a s or az-y- Generically, a s , a^ 7 and a 7 are all present and 
of comparable size, and all three operators can mediate 
four lepton final states. The experimental searches [13, 
l28j require that the invariant mass of the one of the lepton 
pairs is near the Z pole, so the contribution of a 7 is 



and 



azj 



need to be 



small, but we will see that both a s 
considered. 

Bounds on all the operators in Eqs. (P) and @ can be 
set using the absence of single production of this reso- 
nance at LEP [4l[ . In order to interpret these bounds at 
the LHC, however, the production cross section of this 
scalar must be computed, and that is a priori unknown. 

If 4> does couple dominantly via ah that would be ev- 
idence that it is indeed a Higgs. On the other hand, 
it could still be something more exotic such as a dila- 
tion (421 - 1441 ] or a radion [iq . The crucial point is that if 
we are going determine if <p gives mass to the Z boson, 
we must show that its coupling to ZZ* is dominantly 
through ah- As we will show in this paper, the kinematic 
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distributions of the four lepton events can discriminate 
ah from a s and az-y- 

The question of distinguishing these operators via kine- 
matics was considered briefly in (46|. We here extend 
their analysis by studying all possible kinematic variables 
which can distinguish different possible decay operators. 
Furthermore, the analysis of [46[ only considers a mass 
of a scalar greater than twice the Z mass, while we here 
will be working in the kinematic regime where one of the 
Z's is far off-shell. In [3B|, they use kinematic methods to 
distinguish two different kinds of parity even scalars, but 
both of their possibilities are still responsible for giving 
mass to the W and Z. In other words, they both have 
significant a/> in the language of Eq. ([1} . 

The analysis of Q shows that with a fit of the 77, ZZ* , 
and WW* rates, as well as the absence of a large anomaly 
in continuum Z*y, that the scenario of the four lepton 
decays being due to a s is strongly disfavored. As already 
mentioned, while this statement has few assumptions, 
it is still model dependent and we seek to confirm this 
exclusion by a direct measurement. 

If the new resonance has anomalous couplings to Z 
and W , then the production of the resonance through 
vector boson fusion (VBF) would also be modified. These 
effects were studied in (ZJ I47l - l49| . where it was shown 
that angular correlations between the two tagging jets in 
VBF can constrain the value of the operators in Eq. ([T])0 
Here, we only consider decays of the resonance, but these 
two types of measurements can be complementary in fully 
characterizing the nature of the new state. 

The organization of this paper is as follows: in Sec. |TT] 
we describe the kinematics of four lepton events, as well 
as our event generation and selection cuts. We also show 
the distributions which will allow discrimination of the 
different scenarios. In Sec. IIIII we describe the statistical 
procedure we use to distinguish the different scenarios, 
and in Sec. II VI we conclude. 



II. FOUR LEPTON EVENTS 

The kinematics of four lepton events are described in 
detail in many places in the literature, see for exam- 
ple [36|, HH . Here we describe only the variables relevant 
to our analysis. Because we are trying to distinguish dif- 
ferent scalar scenarios, information from the production 
vertex is lost, and angles relative to the beam are irrele- 
vant. The kinematic variables sensitive to decays of the 
scalar are: 

• $ - The angle between the decay planes of the two 
Z bosons in the rest frame of the scalar. 



Monte Carlo code which simulates anomalous VBF events can 
be found in [Holl . 



• 9\ - The angle between the fermion coming from 
the decay of Z\ and the Zi momentum in the Z\ 
rest frame. 

• 62 - 9 1 with Z's interchanged. 

• Mj - The invariant mass of the two Z's. We take 
the convention Mi > M2. 

The distributions for B\ and 62 are the same because of an 
exchange symmetry between the two Z's when they both 
decay to leptons, but it is important to use both because 
that increases the number of observables for each event. 
In most events, M x ~ Mz regardless of the coupling to 
the Z. These variables are all independent subject to 
the constraint (Mi + Ma) < y/s where s is the invariant 
mass squared of the four lepton system. We consider 
four electron, four muon, and two electron + two muon 
events. In the first two types of events there is a two 
fold ambiguity in assigning leptons to parent Z's, which 
we break by taking Mi to be the pair of opposite sign 
leptons whose invariant mass is closest to the Z. 

We compute tree level analytic expressions for the full 
differential decay width in terms of frame invariant 4- 
vector dot products before choosing the frame with kine- 
matic variables described above. These expressions, at 
least of the Higgs, can be found in many places in the lit- 
erature including (38| and references therein. NLO and 
finite mass corrections have been computed in the case of 
the Higgs and have been shown to be a few per cent. 
We plot normalized one dimensional distributions for 
cos 0i, and M2 in Fig.[T] We take = 125 GeV here and 
throughout. Distributions will change little with varia- 
tions of within the experimental resolution. We re- 
strict Mi to the ranges described below in order to better 
approximate the experimental searches. We consider the 
operators c%, a s and az-y turning on one operator at a 
time with the others set to zero. 

In general, all three operators, ah, a s , and az-y will be 
non-zero. In the case of a Higgs-like state which gives 
mass to the Z, a s and az-y are loop suppressed and ah 
will dominate. If the new state does not contribute to 
electroweak symmetry breaking, then ah will often be 
negligible. If a s ~ az-y, which is typically the case if 
the two operators are generated by loops of electroweak 
charged matter, then the effects of a^ 7 will dominate. 
This is because in the allowed region for M2, a photon 
will be much closer to on-shell than a Z, and because 
the Z has suppressed couplings to leptons relative to the 
photon. Even if a s is ten times larger than az-y, azy 
will dominate the decay and we can consider turning on 
just azy as a reasonable approximation. On the other 
hand, one could imagine a model where az-y is very small, 
possibly due to tuning, and we therefore consider turning 
on only a s as a stand in for this possibility. From this 
analysis, we see that in most of the parameter space, one 
operator will dominate over the other two, which is why 
we consider scenarios where only one operator is turned 
on at a time. 
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FIG. 1. Normalized distributions for $ (top), cos#; (middle), 
and M2 (bottom) for m$ — 125 GeV. Each plot shows curves 
from our three different scenarios with ah blue (solid), a 3 red 
(dashed), and az-y green (dot-dashed). 



The top two panels in Fig. Q] show that the angular 
distributions, particularly that of cos 6 provide good dis- 
criminating power between a Higgs-like scenario , and 
the two non-Higgs-like possibilities. The third plot shows 
that the A/ 2 distribution is different for all three scenar- 
ios, and the difference is even more pronounced for small 
values of M^. This can be seen from the following simple 
analysis. For ah, the matrix element goes to a constant as 
Mi — > 0, and a phase space factor of M2CLM2 makes the 
rate go to zero. For a s , the matrix element goes as M% 
for small M.% because of the derivative in the operator, 
so dr falls as M| . Finally, for az-y, the matrix element 
goes as I/M2 because the photon propagator in the de- 
nominator and the derivative in the numerator, and thus 



the rate goes as I/M2. As we will see below, realistic 
detector cuts such as those on lepton px will change this 
low M2 behavior, but this simple analysis shows that if 
the experiments could push down the M 2 reach of the 
events, they would gain discriminatory power. 

We do not include a plot for Mi because in all sce- 
narios, it looks similar with a large peak at Mz that has 
width of Y z- The M\ distribution does, however, provide 
some discrimination power in that the number of events 
well below Mz differs for our three different scenarios. 
For example, in the ah scenario, 70% of the events will 
lie more than 2Y z away from Mz, while the correspond- 
ing fraction for a s {az-y) is 64% (84%). The majority of 
these non-resonant events have Mi < Mz- 

If the four lepton events are dominated by az-y, then 
there should also be decays to on-shell photons. It has 
been pointed out that searching for the Higgs in decays 
to Z7 is a promising channel [53. While there is as yet 
no direct limit in this channel, [3j uses the measurement 
of the cross section to place a limit on the ratio of the 
Z7 mode to the four lepton mode to be about 40. Given 
this, we take the Z"f mode to be an unlikely possibility, 
but we still believe in checking the data to see if it can 
be directly excluded. 

In order to compare to experiment, we also generate 
Monte Carlo (MC) events. We use the Johns Hopkins 
MC described in [351 ] to simulate ah and a s , and Mad- 
graph 5 [EH for az-y- We generate gg — > — > A£ events 
where t = e,n at the LHC with y/s = 8 TeV. Gluon 
fusion is the dominant mode of Higgs production at the 
LHC 
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Since our variables are mostly sensitive to de- 
cay and not production, the errors introduced by ignor- 
ing sub-dominant production modes will be small. We 
require our events to contain four charged leptons (e or 
/i) with 

• p T > 10 GeV 

• |r?| < 2.5 

• 50 GeV < Mi < 110 GeV 

• M 2 > 15 GeV, 

which roughly mimics the experimental selection criteria 
in [13, HH . Histograms for the distinguishing kinematic 
variables from generated events are overlaid on the ana- 
lytic results in Figs. [2] and [3] Because the experimental 
resolution for energy and direction of leptons is so pre- 
cise, we do not apply any smearing to the events. While 
a truly realistic study will need to take into account ex- 
perimental reality, we here see how far the experiments 
could get with just the geometric cuts above. 

In Fig. [21 we plot the cos6*i and cos O2 distributions for 
1000 generated Monte Carlo events which pass the above 
cuts. We compare it to the theoretical distribution which 
is the same for the two angles. We see that the cuts have 
limited effect on cos#i, but the rate for cos 6*2 ~ ±1 is 
suppressed. This is because in that configuration, one 
of the leptons is nearly aligned with the boost direction 
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FIG. 2. Normalized distribution for cosd in the ah scenario. 
The blue (solid) curve is the same as the theory curve from 
Fig-HJ the red (dashed) histogram is the distribution for cos 9i 
for 1000 Monte Carlo events, while the green (dot-dashed) 
histogram is cos 62 for the same events. 
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FIG. 3. Normalized M2 distributions. The blue (solid) curve 
is the theory prediction in the ah scenario, while the light 
blue (dot-dashed) histogram is 1000 Monte Carlo events also 
in the ah scenario. The red (dashed) histogram is 1000 events 
in the a s scenario. 



needed to go to the lab frame from the Z2 rest frame, 
and thus preforming that boost will reduce its energy 
and make it less likely to pass the px cut. This effect is 
small for cos 0\ because the lepton energies in the Z\ rest 
frame are much larger. 

In Fig. |31 comparing the blue (solid) curve to the light- 
blue (dot-dashed) histogram, we see that the experimen- 
tal cuts reduce the event rate for small M^. Even after 
these cuts, however, the histograms for ah and a s still 
differ, so the experimental cuts do not wash out the dis- 
criminating power. 



III. DISTINGUISHING OPERATORS 

In order to estimate the ability of the LHC to discrim- 
inate a Higgs-like scenario dominated by ah from other 
scenarios, we employ a likelihood analysis of the gen- 
erated events. We consider only signal events because 



requiring the invariant mass of the four lepton system to 
be near the mass of the new boson can make the signal to 
background ratio significantly larger than one. Further- 
more, reweighting techniques such as the one laid out 
in [55| can be used to further purify the event selection. 

We use a standard unbinned likelihood analysis which 
is described in detail in 35]. We can use the computed 
normalized differential cross section as a probability dis- 
tribution P($, 9i, Mi\ai) for each operator ah, a s , and 
az-y- The normalization is computed with the Mi cuts 
described above because they are independent of Lorentz 
frame. Taking the px and 77 acceptance into account in 
P would improve the statistical power of the test, but 
because those cuts are frame-dependent, we leave that to 
further work. 

Given a sample of N events, we can then construct a 
likelihood C(di) — n^Li Pj( a i)- With this likelihood we 
can then compare two different scenarios, a\ and ai by 
constructing a hypothesis test with test statistic defined 

by m 



A = 21og[£(oi)/£(a 2 )]. 



(3) 



Since we are taking the resonance mass as input and us- 
ing the normalized differential cross sections to construct 
our likelihood functions, there are no free parameters 
(nuisance parameters) in this ratio, making this a simple 
hypothesis test. 

To estimate the expected significance of discriminating 
between two different hypotheses corresponding to two 
different operators, we follow a similar analysis to that 
found in [351 ] . To begin, we take one hypothesis as true, 
say ai and generate a fixed number N of a\ events. We 
then construct A as above for a large number of pseudo- 
experiments each containing N events in order to obtain 
a distribution for A. We then repeat this exercise tak- 
ing ai to be true and again obtain a distribution for A. 
These two distributions are shown in Fig. [4] comparing 
ah and a s . This figure shows 5000 pseudo-experiments of 
50 events each, which shows a clear separation between 
the two scenarios. 

With the two distributions for A in hand we can com- 
pute an approximate significance by the following proce- 
dure. If we denote the distribution with negative mean 
as / and the distribution with positive mean as g, we find 
a value A such that 



fdx = 



gdx. 



(4) 



Schematically, this value of A corresponds to a value such 
that if the experiment observed that value for the test 
statistic, it would have no discriminatory power between 
the two scenarios. We then interpret the probability 
given by either side of Eq. (|1]) as a one sided Gaussian 
probability, which can then be interpreted in terms of 
number of a. This procedure is shown schematically in 
Fig. 2] with the areas of the two shaded regions being 
equal and corresponding to the probability of excluding 
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FIG. 4. Normalized distribution of our test statistic A when 
a/, is true on the right (blue), and when a s is true on the 
left (pink). Each histogram is the result of 5000 pseudo- 
experiments with 50 events each. The vertical (green) line 
is A defined in Eq. (j4]) such that the area to the right of A 
under the a„ histogram is equal to the area to the left of A 
under the histogram. We also draw a Gaussian over each 
histogram with the same median and standard deviation. 



the correct hypothesis. For a simple hypothesis test, this 
Gaussian approximation is often sufficient [56| . and we 
see from Fig. [4] that the A distributions are well approx- 
imated by Gaussians. 

This procedure is repeated many times for a range of 
numbers of events N to obtain a significance as a func- 
tion of N for each hypothesis. We show this for the case 
where a\ = ah and a-i — a s or a-i — az- f in Fig. [5] We 
see that with O(50) events, we can distinguish renormal- 
izable from nonrenormalizable coupling to ZZ at 95% 
confidence, and with 0(100) events we can get a 99% ex- 
clusion. The operator azj can be distinguished from ah 
at 95% confidence with as few as 20 events. The third 
possibility, which we do not show, is even easier; a s and 
az-f can be distinguished from one another at 95% with 
just 10 events. 



IV. CONCLUSIONS AND OUTLOOK 

Testing the properties of the newly discovered reso- 
nance near 125 GeV is of utmost importance. While the 
rate and branching ratio data are consistent with the new 
particle being the Standard Model Higgs, direct tests of 
its properties are still essential. In this paper we have 
examined the discriminating power of events where the 
new particle decays to four leptons. These events can be 
used to measure the Lorentz transformation properties 
of this particle, but even if it is confirmed to be a par- 
ity even scalar, it still need not be the Higgs; it could 
couple to the gauge bosons via higher dimensional oper- 
ators rather than via the renormalizable operator in the 
Standard Model. 

We have analyzed how well kinematic distributions in 
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FIG. 5. Expected significance as a function of number of 
events in the case of a^ vs a 3 on top, and vs az~, on bot- 
tom. We use a different horizontal scale for the top and bot- 
tom plots because far fewer events are needed to discriminate 
from az-y than from a B . We also fit with a function pro- 
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four lepton events can distinguish between different ten- 
sor structures of the coupling to gauge bosons. In par- 
ticular, we looked a coupling directly to Z^Z^ 1 , as well 
as couplings to a pair of field strength tensors of the 
Z, and a coupling to the field strength of the Z and 
of the photon. All three scenarios will produce one lep- 
ton pair near the Z pole, while the other pair will have 
much lower invariant mass. We find that with O(50) sig- 
nal events, a Higgs-like state can be discriminated from 
ZZ field strength tensor couplings with 95% confidence, 
while only 20 events are needed to make the same deter- 
mination for field strength coupling to Zj. This shows 
that the 2012 LHC run has excellent prospects to con- 
strain the tensor structure of the new state's coupling to 
gauge bosons. 

While the four lepton final state is one of the most 
powerful for discriminating different scenarios, it would 
be interesting to look at kinematic variables in other final 
states. For example, in the decay to WW* where both 
W's decay leptonically, the angles between the leptons 
and the transverse angles with missing energy will pro- 
vide discriminating power, though this channel is difficult 
because of the large background. A search for decay to 
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where the photon is on-shell would give a direct mea- 
surement of the az~f coupling to the Zj field strength 
operator given in Eq. ([2]l. If this mode is in fact ob- 
served, kinematic analysis of final states in that channel 
could further uncover the nature of the new particle. 

There are many ways to both directly and indirectly 
learn about the couplings of the new state. We have ar- 
gued that even if the state is a parity even scalar, it could 
be decaying to four leptons in a very non-Higgs-like way, 
possibly even through Z7* instead of ZZ*. Therefore, 
we hope that this work will inspire new measurements 
which can have strong discriminating power in the very 
near future. 



Note Added: While completing this paper, we received 
Refs. [57J, |58| which discuss similar ideas. 

We thank Zackaria Chacko, Kyle Cranmer, Andrei Grit- 
san, Ian Low, and Kirill Melnikov for useful discussions. 
We also thank Markus Schulze and Nhan V. Tran for 
their help with the Johns Hopkins Monte Carlo, and 
Kunal Kumar for help with Madgraph. DS is sup- 
ported in part by the NSF under grant PHY-0910467 
and gratefully acknowledges support from the Maryland 
Center for Fundamental Physics. DS is also supported 
in part by the Project of Knowledge Innovation Pro- 
gram (PKIP) of Chinese Academy of Sciences, Grant No. 
KJCX2.YW.W10. 



S. Chatrchyan et 



[CMS Collaboration], 



[arXiv: 1207.723 5 [hep-ex]]. [27 

G. Aad et al. [ATLAS Collaboration], [arXiv:1207.72T4l 

[hep-ex]]. [28; 

I. Low, J. Lykken and G. Shaughnessy, arXiv: 1207. 1093 
[hep-ph]. [29 
T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile and 
M. C. Gonzalez-Garcia. larXiv:1207.1344l [hep-ph], [30 
A. Arbey, M. Battaglia, A. Djouadi and F. Mahmoudi, 
larXiv:1207.1348l [hep-ph]. 

P. P. Giardino, K. Kannike, M. Raidal and A. Strumia, [31 

larXiv:1207.1347l [hep-ph], 

M. R. Buckley and D. Hooper. larXiv: 1207. 14451 [hep-ph] . 
J. Baglio, A . Djouadi and R. M. Godbole, [32 

larXiv: 1207. 14511 [hep-ph], 

J. F. Gunion, Y. Jiang and S. Kraml, larXiv:1207.1545l [33 
[hep-ph]. 

J. Ellis and T. You JarXiv:f2 07.1693 [hep-ph]. [34 
M. Montull and F. Riva. larXiv:1207.1716l [hep-ph], 
J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott, [35 
larXiv: 1207. 17171 [hep-ph]. 

D. Carmi, A. Falkowski, E. Kuflik, T. Volansky and 

J. Zupan, larXiv: 1207.T7181 [he p-ph], [36 

A. G. Cohen and M. Schmaltz, larXiv: 1207.34951 [hep-ph] . 
S. Banerjee, S. Mukhopadhyay and B. Mukhopadhyaya, 

larXiv:1207.3588l [hep-ph], [37 

J. Cao, Z. Heng, J. M. Yang and J. Zhu JarXiv:1207.3698l 

[hep-ph], [38 

D. Bertolini and M. McCullough, larXiv: 1207.42 09 [hep- 
ph]. [39 
A. Joglekar, P. Schwaller and C. E. M. Wagner, 
larXiv: 1207.42351 [hep-ph], [40 
N. Arkani-Hamed, K. Blum, R. T. D'Agnolo and J. Fan, 
larXiv: 1207.44821 [hep-ph] . [4 1 

F. Bonnet, T. Ota, M. Rauch and W. Winter, 

larXiv:1207.4599l [hep-ph], [42 

N. Craig and S. Thomas, larXiv: 1207.48351 [hep-ph] . 
S. Matsuzaki and K. Yamawaki, larXiv:1207.5911l [hep- [43 
ph]. 

J. Moffat, larXiv:12 07.6015 [hep- ph], [44 
T. Plehn and M. Rauch, larXiv: 1207.61081 [hep-ph] . 
I. Dorsner, S. Fajfer, A. Greljo and J. F. Kamenik, [45 

larXiv:1208.1266l [hep-ph], 

[26] B. Coleppa, K. Kumar and H. E. Logan. larXiv:1208.2692l [46 



[hep-ph]. 

CMS Collaboration, CERN, July 2012, 

http : / /cdsweb . cern . ch/ record/ 1460664 

ATLAS Collaboration, ATLAS-CONF-2012-092, CERN, 

July 2012 http://cdsweb.cern.ch/record/1460411 

A. Soni and R. M. Xu, Phys. Rev. D 48, 5259 (1993) 

[hep-ph/9301225] . 

-m. Cheung, A. Djouadi, B. A. Kniehl 
Phys. Rev. D 49, 79 (1994) 



V. D. Barger, K. 
and P. M. Zerwas, 
[hep-ph/9306270] . 
S. Y. Choi, D. J. 
and P. M. Zerwas, 



Miller, 2, M. M. Muhlleitner 
Phys. Lett. B 553, 61 (2003) 
[hep-ph/02 10077] , 

C. P. Buszello, I. Fleck, P. Marquard and J. J. van der 
Bij, Eur. Phys. J. C 32, 209 (2004) [hep-ph/02 12396] , 
W. -Y. Keung, I. Low and J. Sh u, Phys. Rev. Lett. 101, 
091802 (2008) [arXiv:0806. 28641 [hep-ph]]. 

0. Antipin and A. Soni, JHEP 0810, 018 (2008) 
|arXiv:0806.3427l [hep-ph]] . 

Y. Gao, A. V. Gritsan, Z. Guo, K. Melnikov, M. Schulze 
and N. V. Tra n, Phys. Rev. D 81, 075022 (2010) 
[arXiv:1001. 33961 [hep-ph]]. 

A. De Rujula, J. Lykken, M. Pierini, C. Rogan 
and M. Sp iropulu, Phys. Rev. D 82, 013003 (2010) 
|arXiv:1001.5"300l [hep-ph] ] . 

U. De Sanctis, M. Fabbrichesi and A. Tonero, Phys. Rev. 
D 84, 015013 (2011) |arXiv: 1103.19731 [hep-ph]]. 
J. S. Gainer, K. Kumar, I. Low and R. Vega-Morales, 
JHEP 1111, 027 (2011) 1 arXiv: 1108.22741 [hep-ph]]. 

1. Low and J . Lykken, JHEP 1010, 053 (2010) 
|arXiv:1005. 08721 [hep-ph]] . 

P. J. Fox, D. Tucker-Smith and N. Weiner, JHEP 1106, 

127 (2011) |arXiv:1104.5450l [hep-ph]]. 

V. Hankele, G. Klamke, D. Zeppenfeld and T. Figy, Phys. 

Rev. D 74, 095001 (2006) |hep-ph/0609"075] . 

E. Gildener and S. Weinberg, Phys. Rev. D 13, 3333 

(1976). 

W. D. Goldberger, B. Grinstein and W. Skiba, Phys. Rev. 
Lett. 100, 111802 (2008) |arXiv:0708. 14631 [hep-ph]] . 
J. Fan, W. D. Goldberger, A. Ross and W. Skiba, Phys. 
Rev. D 79, 035017 (2009) |arXiv:0803.2040l [hep-ph]]. 
W. D. Gold berger and M. B. Wise, Phys. Rev. Lett. 83, 
4922 (1999) [hep-ph/9907447] . 

Q. -H. Cao, C. B. Jackson, W. -Y. Keung, I. Low and 



7 



J. Shu, Phys. Rev. D 81, 015010 (2010) |arXiv:0911.3398l 
[hep-ph]]. 

[47] T. Plehn, D. L. Rainwater and D. Zeppenfeld, Phys. Rev. 

Lett. 88, 051801 (2002) [hep-ph/0105325] . 
[48] G. Klamke and D. Zeppenfeld, JHEP 0704, 052 (2007) 

[hep-ph/0703202] [HEP-PH] ] . 
[49] K. Hag iwaraTOT Li and K. Mawatari, JHEP 0907, 101 

(2009) [arXiv:0905.4"3l4l [hep-ph]] . 
[50] K. Arnold, J. Bellm, G. Bozzi, F. Campanario, 

C. Englert, B. Feigl, J. Frank and T. Figy et ai, 

larXiv:1207.4975l [hep-ph], 

[51] B. A. Kniehl and O. L. Veretin.i arXiv:1206.7110l [hep-ph]. 
[52] J. S. Gainer, W . -Y. Keung, I. Low and P. Schwaller, 

larXiv:1112.1405l [hep-ph]. 
[53] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and 

T. Stelzer, JHEP 1106, 128 (2011) |arXiv: 1106.05221 



[hep-ph]]. 

[54] S. Dittmaier et al. [LHC Higgs Cross Section 
Working Group Collaboration], I arXiv: 1101. 05931 [hep- 
ph]. See also https://twiki.cern.ch/twiki/bin/view/ 
LHCPhysics/CrossSections 

[55] M. Pivk and F. R. Le Diberder, Nucl. I nstrum. Meth. A 
555 (2005) 356 |arXiv:physics/0402083| . 

[56] R. Cousins, J. Mumford, J. Tucker and V. Valuev, JHEP 
0511, 046 (2005). 

[57] S. Bolognesi, Y. Gao, A. V. Gritsan, K. Mel- 
nikov, M. S chulze, N. V. Tran and A. Whitbeck, 
larXiv:1208.40l8l [hep-ph]. 

[58] R. Boughezal, T. J. LeCompte and F. Petriello, 
larXiv:1208.43TT1 [hep-ph] . 



